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I. INTRODUCTION 


If one were to ask a randomly chosen physicist 
oscillator 


whether the equation of motion of the damped classical harmonic 


X -b 27X -b W^X = 0 (7 > 0) 


( 1 ) 


can be derived from a time-independent Hamiltonian, almost certainly the answer would be a resounding “no,” because 
multiplying this equation by x, one obtains the equation 


dt 




^CU X 


—27 x ^. 


( 2 ) 


The quantity ^x^ -b jw^x^, which appears to be the sum of a kinetic and a potential energy, is not conserved and 
decreases with time. So, one might think that © cannot be derived from a time-independent Hamiltonian. 

However, Bateman yJl made the remarkable observation that if one appends the time-reversed oscillator equation 
with undamping (gain) instead of damping, 

y - 2jy-{-uj'^y = 0 (7 > 0), (3) 


then even though the two oscillators are independent and noninteracting, the two equations of motion o and ([3]) can 
be derived from the time-independent quadratic Hamiltonian 

H =pq + j{yq - xp) -b - 7^) xy. (4) 


The two oscillator equations follow directly from Hamilton’s equations of motion 


dH 


dp 
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( 5 ) 

( 6 ) 

( 7 ) 

( 8 ) 


To derive © we differentiate 
to derive @, we differentiate 


([5|) with respect to t, eliminate q by using ([5]), and eliminate q by using 
with respect to t, eliminate p by using ©, and eliminate p by using 


®. 

®. 


Similarly, 
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The Hamiltonian 0 is VT symmetric 0; under parity reflection V the oscillators with loss and gain are inter¬ 
changed, 

V : X ^ y, y ^ X, q, q ^ p, (9) 

and under time reversal T the signs of the momenta are reversed, 

T-.x^x, y^y, p^-p, q^-q- (10) 

[The Hamiltonian i? in (|3]) is VT symmetric but it is not invariant under V or T separately.] Because the balanced- 
loss-gain system is described by a time-independent Hamiltonian, the energy (the value of H) is conserved in time. 
However, the energy has the complicated form in (|3]) and is not a simple sum of kinetic and potential energies. If the 
loss and gain terms in o and ([3]) were not exactly balanced (that is, if the velocity terms did not have equal but 
opposite signs), the system would not be derivable from a time-independent quadratic Hamiltonian. 

It is even more remarkable that the equation of motion o of the damped oscillator can be derived from a (non¬ 
quadratic) time-independent Hamiltonian without havirw to introduce any additional degrees of freedom. The con¬ 
struction of such a Hamiltonian was first given in Ref. (3| and was accomplished by using a rather obscure technique 
called the Prelle-Singer method. 

In this paper we show how to construct the Hamiltonian for an arbitrary homogeneous linear constant-coefficient 
differential equation of any order. First, in Sec. [TT] we do so for the second-order equation ([T]) and in Sec. [ITT] we 
demonstrate the procedure for a general third-order equation. An interesting special case of such an equation is 
the equation that describes the nonrelativistic self-acceleration of a charged oscillating particle and it is quite 
remarkable that even though there are runaway modes, the energy of such a system is conserved. Then, in Sec. IIVI 
we generalize our procedure to an arbitrary nth-order constant-coefficient equation. Section |V] discusses the problem 
of quantization and we show that quantizing the classical Hamiltonians discussed in this paper is nontrivial. Finally, 
Sec. EH gives a brief summary. 


II. HAMILTONIAN FOR A GENERAL LINEAR CONSTANT-COEFFICIENT SECOND-ORDER 

DIFFERENTIAL EQUATION 


We begin with m and substitute x{t) = e 


This gives a quadratic equation for the frequency n: 


-I- 2iji' — = 0. 


This equation factors 


(n - wi) {v - 0J2) = 0 , 


where 


wi -I- a;2 = —2^7, wia;2 = — 


and thus 


101^2 = —*7 ± H = —i'y ± 




— 7^. 


(II) 


( 12 ) 


(13) 


(14) 


We claim that a Hamiltonian H that generates a general linear constant-coefficient nth-order evolution equation 
has the generic form 


H = axp + f{p), 


(15) 


where a is a constant and f{p) is a function of p only. For the case of the second-order equation ([T]), one such 
Hamiltonian is 

Hi = -iujixp + ^ (16) 

UJi — UJ2 

A second and equally effective Hamiltonian is obtained by interchanging the subscripts 1 and 2: 


H 2 = —iuj2xp + 




a;2 — wi 


(17) 





3 


These Hamiltonians appear in Ref. Q for the case of over-damping ( 7 ^ > in which case they are real, but they 
apply equally well when ( 7 ^ < (We are not concerned here with the reality of the Hamiltonian.) 

For the Hamiltonian Hi, Hamilton’s equations read 


dHi , 

X = —— = —iuiix+p \ 
op 


p = 


dHi 

dx 


= lOJip. 


(18) 

(19) 


We then take a time derivative of mi) and simplify the resulting equation first by using mi) and then by using m-- 

X + iujix = — — 


Wi 


— ix}2ltOi 


Thus, 


= —IW2P 
= —iuj 2 (x + iuJix). 


X + i (oJi -b 012 ) X — UJ1UJ2X = 0, 


( 20 ) 

( 21 ) 


which reduces to m upon using (US- 

The evolution equation © has one conserved (time-independent) quantity, and this quantity can be expressed in 
terms of the function x(t) only. To find this quantity, we begin with (1181) and solve for p: 

p = {x + iuJix)~‘^^^‘^^ . ( 22 ) 

We then use this result to eliminate p from the Hamiltonian Hi. Since Hi is time-independent, we conclude that 


(i + ioj2xY'^ 
(x + iuJix)‘^^ 


(23) 


is conserved. Had we started with the Hamiltonian H 2 we would have obtained the conserved quantity 


_ (i -b iuJixY^ 

^ {x + iuj 2 x)‘^^ ’ 


(24) 


but this is not an independent conserved quantity because C 2 = l/Ci. These conserved quantities were also found in 
Ref. 0 for the case of over-damping. 

When 7 = 0 , these results reduce to the familiar expressions in the case of the simple harmonic oscillator. In this 
case we let oj = cui = —uj2 so that Hi becomes 


Hi = -iu}xp+ (25) 

which is related to the standard simple harmonic oscillator Hamiltonian by the change of variable p —>■ p — iuix. The 
conserved quantities C 2 and Ci become simply (i^ -ba;^a;^)^“, in which we recognize the usual conserved total energy. 


III. HAMILTONIAN FOR A CONSTANT-COEFFICIENT THIRD-ORDER EQUATION 

In this section we show how to construct a Hamiltonian that gives rise to the general third-order constant-coefficient 
evolution equation 


{D + iuJi){D + iuj2){D + ioj3)x = 0, (26) 

where D = The Hamiltonian that we will construct has just one degree of freedom. 

An interesting physical example of such a differential equation is the third-order differential equation 

mx + kx — mr'x = 0 (27) 

that describes an oscillating charged particle subject to a radiative back-reaction force Q. Following Bateman’s 
approach for the damped harmonic oscillator, Englert Q showed that the pair of noninteracting equations dUl) and 

my + ky -\- mr y = 0 (28) 
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can be derived from the quadratic Hamiltonian 


ps — rq 2rs 
H = - - - H-^ 


pz + qw 


rriT 


rriT'^ 


mzw 

— -h kxy. 


(29) 


This Hamiltonian contains the four degrees of freedom {x,p), {y,q), (z;,r), and {w,s). An interacting version of this 
model was studied in Ref. In fact, we find that the two equations of motion (l?71) and (E51) can be derived from 
the simpler quadratic Hamiltonian 


H = 


pr + qz 



-h kxy^ 

T 


which has only the three degrees of freedom {x,p), {y,q), and {z,r). A similar three-degree of freedom Hamiltonian 
was also found in Ref. 

Our objective here is to find a one-degree-of-freedom Hamiltonian that can be used to derive the third-order 
differential equation (l26l) . Note that the general solution to (l26l) is 

X = -b 026 -“^* -b 036 -*“^^*, (30) 


where ak are arbitrary constants. If we form {D + iuj 2 ){,D -b iuJ^)x, that is, x + i{uJ 2 -b 0 J 3 )x — u! 2 i^ 3 X, we obtain 

= jP + + iuj3)x 

(wi - a;2)(a;i - wa) 

in which the constants 02 and 03 do not appear. Similarly, we have 

{D + iuJ3)(D -b iu}i)x 


026 




a^e 




(012 - W3)(a;2 - wi) 

[D -b iuJi){D -b iuj2)x 
{UJ3 - U}i){uJ3 - UI2) 


(31) 


(32) 


So, assuming that the frequencies ujk are all distinct, there are two independent conserved quantities, namely 

^ _ [x + i(wi -\-W2)x — UJiIjJ2xY^‘^^ 

^ [x + i{ijJ 2 + UJ3)x — 

^ _ [x + i{u>i + i03)x — WiU}3xY/‘^'^ 

^ [x + i{oJ 2 + 0J3)x — 

These expressions and the equation of motion can be derived from the Hamiltonian 

TJ ■ I ^2^03 I ^ 3^1 1 —CJs/cJi 

H = —ILOiXp-\ - p ^-I- p \ 

Wl — UJ2 UJl — UI3 

where 62 and 63 are arbitrary constants. Thus, p = — ^ = iuJiP- This means that p oc so that 1/p is directly 

related to the combination in (IHTl) . 

Then, from Hamilton’s equation x = ^ and from further differentiation with respect to t, we obtain 

X = -iwia;-b -b 

X = -iuJix - iuj2b2P~‘^^^‘^^ - iuj3b3P~‘^’^^‘^\ (35) 

X = —iwii — 


(33) 

(34) 


These equations depend on the constants &2 and 63. Nevertheless, after we combine these equations and perform some 
simplifying algebra, we obtain 


'x + iigUi -b a ;2 -b UJ 3 )x — {uJiLU2 -b Ul2<-03 + 0 J 3 UJi)x — i 0 Ji(jJ 20 J 3 X = 0. (36) 

The constants &2 and 63 have disappeared in this combination and we have reconstructed the equation of motion (1261) . 
Using only derivatives up to the second order, we can find expressions for 62 ^““’^^“’^ and b 3 p~'^^/^^, namely 

i{uj3 — UJ2)b2P~‘^'^^‘^^ = X + i{uJi + UJ3)x — OJ1UJ3X, 
i{uj2 — UJ3)b3P~‘^^^‘^^ = X + i{u!i + UJ2)x — LO1UJ2X. 


( 37 ) 
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These are precisely the combinations appearing in p3l) . and from them we can construct the conserved quantity 
C2/C3. 

In order to derive the second conserved quantity we use the fact that the Hamiltonian is a constant. We then 
evaluate H in terms of a;, x, and x using (1371) and after some algebra we find that 


H = iujip 


X + i{(jJ2 + UJ3)x — UJ2OJ3X 
(wi - UJ2)iuJl - UJ3) 


(38) 


which contains precisely the other linear combination of x, x, and x that appeared in m- As already mentioned, 
1 Ip is proportional to that combination. 

To summarize, the evolution equation (I26|) can be derived from the unusual time-independent Hamiltonian (1341) 
containing the single coordinate variable x and its conjugate momentum p. This Hamiltonian is a conserved quantity, 
which can be expressed as a function of x, x, and x. There is a second conserved quantity having a similar structure. 
In the next section we show how this procedure generalizes to a linear constant-coefficient differential equation of 
order n, for which there are n — 1 conserved quantities involving time derivatives up to order n — 1. 

However, before moving on, it is important to resolve an apparent paradox, namely, how a Hamiltonian with a 
single degree of freedom can give rise to a differential equation whose order is greater than two. The problem is this: 
Our Hamiltonian has the general form in (1151) : H = axp + f{p). Therefore, the equations of motion are simply 


X = ax -I- gip), (39) 

where g{p) = f'{p), and 

p = —ap (40) 

We solve (|40)) first, 

p{t) = Ce-“‘, (41) 

where C is an arbitrary constant. Next, we return to ([39]), which becomes 

X = ax + g (42) 

after we eliminate p by using (TO- This is a first-order equation. Thus, its solution has only two arbitrary constants: 

x(t) = (j}{t,C, D). (43) 

We obtained the higher-order differential equation (l30l) by the sequence of differentiations in (1361) . However, the 
solution to an nth-order equation can incorporate n pieces of data such as n initial conditions: x(0), x(0), x(0), x (0), 
and so on. How is it possible to incorporate n pieces of data with only two arbitrary constants C and D7 There 
appear to be n — 2 missing arbitrary constants. 

The surprising answer is the n — 2 pieces of initial data determine n — 2 parameters multiplying each of the fractional 
powers of p in H. (One parameter can always be removed by a scaling.) We call these parameters b^. Thus, we have 
incorporated the initial data into the Hamiltonian in the form of coupling-constant parameters. 

For the triple-dot equation, we can see from (1571) that the ratio is related to the initial conditions. So, 

for the case of the third-order equation, the three arbitrary constants are C, D, and . We emphasize that 

if one just wants a Hamiltonian that gives the equations of motion, the coefficients bk, which play the role of coupling 
constants, are irrelevant and we have shown that they drop out from the equation of motion. However, the coupling 
constants in the Hamiltonian are required to incorporate the initial data and are determined by the initial data. 


IV. HAMILTONIAN FOR A CONSTANT-COEFFICIENT nTH-ORDER EQUATION 


It is straightforward to generalize to the case of an arbitrary nth-order constant-coefficient evolution equation 


whose general solution is 


n 

(U -b iiOr) 

r—1 


x{t) = 0, 


x{t) =''^are 

r—1 


(44) 


( 45 ) 
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For simplicity, we assume first that the frequencies uJr are all distinct; at the end of this section we explain what 
happens if some of the frequencies are degenerate. 

Corresponding to (l32|) and (1^ . we have 


Thus, the quantity 


oc 


(D + iOJr) 

r^s 


x{t). 


Qs = 


(£) + iuJr) 

r^s 


x{t) 




(46) 


(47) 


is proportional to for all s. Hence, the n — 1 independent ratios Qs/Qi (s > 1) are all conserved. Any other 
conserved quantities can be expressed in terms of these ratios. 

The equation of motion and the conserved quantities can be derived from the Hamiltonian 


n 

H = —iujixp + 


UJi — UJr 


(48) 


which is the nth order generalization of (1341) for the cubic case. In this expression the n — 1 coefficients br are arbitrary. 
Note that in constructing the Hamiltonian H there is nothing special about the subscript “1” and it may be replaced 
by the subscript “s” (1 < s <n). 

Degenerate frequencies 

Until now, we have assumed that the frequencies uir are all distinct. However, if some of the frequencies are 
degenerate, there is a simple way to construct the appropriate Hamiltonian: If the frequencies oji and 0 J 2 are equal, 
we make the replacement 

pl-u,2/u:^ 

UJx — LU2 

(In making this replacement we are shifting the Hamiltonian by an infinite constant.) Thus, for uii = UI 2 the 
Hamiltonian Hi va dT51) reduces to 


Hi = —iojixp + logp. (50) 

Hamilton’s equations for this Hamiltonian immediately simplify to (EB) with oJi = (jJ 2 - 
Similarly, for the case wi = UJ 2 the Hamiltonian (13411 reduces to 

H = —iujixp + 62 logp H-^ 52 ^^ 

wi - wa 

and Hamilton’s equations for this Hamiltonian readily simplify to (13611 with uji = UJ 2 - 

Also, if the frequencies are triply degenerate, = W 2 = wa = w, the Hamiltonian in (IMl) is replaced by 

H = —iujxp-\-b\ogp-\- \c{\ogpY, (52) 

where b and c are two parameters that are determined by the initial data. Once again, Hamilton’s equations for this 
Hamiltonian combine to give (l36ll with uti = 102 = 003 = w. 


V. QUANTIZATION 

The obvious question to be addressed next is whether it is possible to use the Hamiltonians that we have constructed 
to quantize classical systems that obey linear constant-coefficient evolution equations. Let us begin by discussing the 
simple case of the quantum harmonic oscillator (QHO), whose Hamiltonian Hi is given in (1^51) . 

One possibility is to quantize Hamiltonian in p-space by setting x = id/dp. Then the time-independent Schrddinger 
eigenvalue equation is, by shifting E hy uj 

HMp) = (^wp^ + ip2^ 


ijj{p) = Eijj{p), 


(53) 
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whose solution is 

^(p) (54) 

In this way of doing things we can derive the quantization condition by demanding that "0 be a well-defined, nonsingular 
function, which requires that E = nuj, where n is a nonnegative integer Q. However, these “momentum-space” 
eigenfunctions are problematical because p has no clear physical interpretation as a momentum, and it is not a 
Hermitian operator. The momentum eigenfunctions are certainly not orthonormal in any simple sense because they 
do not solve a Sturm-Liouville boundary-value problem 

However, we can calculate the corresponding x-space eigenfunctions by Fourier transform using the formula Q 

H„{z) = / dp EP^p^e-P . (55) 

^ V ^ J — CO 


We find that 

i>nix) IX e~2‘^ “ Pn{x), 


(56) 


where (pn{x) is the nth eigenfunction of the QHO. This is consistent with our remark above that Hi is related to the 
standard QHO Hamiltonian by the transformation p ^ p — iujx. This transformation is achieved at the operator level 
by the similarity transformation p —>■ [l^. Because of this additional factor, our eigenfunctions are 

orthonormal with respect to the metric rj = ^ . As an alternative approach, we can cast (051) in x-space as 


Hi 



(57) 


from which we can obtain the ipn (x) directly. 

To summarize, the quantized version of (1251) corresponds to a transformed version of the QHO, where the x-space 
eigenfunctions are simply related to the standard eigenfunctions, and are orthonormal with respect to an additional 
weight function. The p-space eigenfunctions can be written down but their interpretation is not at all obvious (the 
operator p corresponds to the conventional raising operator a^) and are not orthogonal in any simple way. In p space 
the weight function e‘^ “ becomes the highly nonlocal operator a /ap ^ 

If we now generalize to the damped harmonic oscillator, we can still find a solution 0(p) to the time-independent 
Schrodinger equation, namely 


'0(p) oc exp 


Wi 


(Wl - UJ2)‘^ 


l—UJllw2 


(58) 


but even if we take E = noji in order to make the prefactor nonsingular, we are still left with a nonintegral, and 
in general complex, power of p in the exponential. (See, also, the comments in Ref. Q.) Thus, in addition to the 
previously discussed problems with '0(p), we would now have to consider it to be a function in a cut plane. Moreover, 
there is no simple formula like (1551) whereby one could obtain the x-space eigenfunctions. Furthermore, if we cast the 
equation in x-space we obtain 


Hi 


1 

1 — ljJ2/u)i 






(59) 


in which the difficulty associated with a fractional derivative is manifest. 

We conclude that quantizing Hamiltonians of the form in (1151) is nontrivial. The problem of quantizing the cubic 
equation describing the back-reaction force on a charged particle was solved in Ref. (6|. However, the system that was 
actually quantized was a pair of coupled cubic equations in the unbroken T^T-symmetric region. Thus, it may be that 
the most effective approach for quantizing a Hamiltonian of the form (1151) is to introduce a large number of additional 
degrees of freedom. 


VI. SUMMARY 

We have shown that any nth-order constant-coefficient evolution equation can be derived from a simple Hamiltonian 
of the form dUD. Remarkably, this Hamiltonian only has one degree of freedom, that is, one pair of dynamical variables 









{x,p). Furthermore, we have shown that for such a system there are n — 1 independent constants of the motion and 
we have constructed these conserved quantities in terms of x{t). However, we find that it is not easy to formulate a 
general procedure to quantize the system described by the Hamiltonian, and this remains an interesting open problem. 
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